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Abstract 

An impossibility theorem on approximately simulating quantum non-integrable 
Hamiltonian systems is presented here. This result shows that there is a trade- 
off between the unitary property and the energy expectation conservation law in 
time-descretization of quantum non-integrable systems, whose classical counter- 
Sh I part is Ge-Marsden's impossibility result about simulating classically non-integrable 

Hamiltonian systems using integration schemes preserving symplectic (Lie-Poisson) 
O" 1 ' property. 
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^ 1 Introduction 

. P.. 

Recently, much attentions are directed to investigate the interrelation between physics 
and computation. To connect physics with computation, we can classify the problems 
into the the following classes: 



class (1) 

class (2) 

class (3) 

class (4) 



Connection between classical physics and classical computation, 
Connection between quantum physics and classical computation, 
Connection between classical physics and quantum computation, and 
Connection between quantum physics and quantum computation. 



Concerning the class (4), simulating quantum behavior such as quantum chaos using 
classical computers is known to be a notoriously difficult computational problem B. One 
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of the main difficulties is that one must discretize a continuous time parameter of equations 
of motion in order to simulate on computers. Thus, it is an important question to ask 
whether we can always have a suitable time-discretization scheme of Schrodinger equation. 

In this paper, I will give a somewhat negative answer to this question: In the case 
of quantum non-integrable systems with an explicit time-independent Hamiltonian oper- 
ator, no explicit time-discretization algorithm preserving unitary property can simulate 
quantum non-integrable behavior without violating the conservation law of energy expec- 
tation. Since the original quantum nature must have these two properties, namely, the 
conservation law of energy expectation and the unitary property of time evolution, this 
means that there is a fundamental limit in simulating quantum non-integrable behavior 
using unitary maps like quantum computers. This negative result can be regarded as a 
quantum analogue of Ge-Marsden's theorem0: No symplectic integrator can simulate 
non-integrable behavior in a class of autonomous Hamiltonian systems without violating 
the energy conservation. These fundamental limits, whether quantum or classical, suggest 
the importance of the notion of integrability in simulating physical behavior. In Section 
2, we give a brief explanation of time- discretization preserving unitary property. In Sec- 
tion 3, we give a theorem about the impossibility of simulating quantum non-integrable 
systems. In Section 4, we discuss various aspects about our results. 



2 Simulation technology preserving unitary prop- 
erty 

Time-evolution of quantum computation can be seen as a class of successive iterations of 
unitary transformations |l|, |^, ||. The time-evolution operator has the form of 

U(At) = exp [-iAtH/h] , (1) 

where if is a Hamiltonian operator with Hermitian property, At is the time duration of 
each computation process and an exponential operator exp [xA] is defined as 

oo ( T A\n 

exp [xA] = ]T x = -iAt/h. (2) 

n=0 n - 

Let A and B be Hermitian operators as the generators of two different elementary pro- 
cesses of unitary dynamics. In general, A does not commute with B: 

[A, B} = AB - BA^ 0. (3) 

To track computational processes retaining unitary property, evaluating the following 
time-evolution operator 

exp [x(A + B)} (4) 
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is relevant to various problems. In fact, there are infinitely many methods to get pertur- 
bation series of Eq.(^j). The Feynman path-integral method|| 



exp [x(A + B)\ 



1 + 



x(A + B) 



n 



(5) 



discovered in his study of quantum electro-dynamics is a first-order method based on the 
identity 

x(A + B) 



exp [x(A + B)\ = lim 



1 + 



n 



(6) 



However, the above approximation breaks unitary property in each elementary dy- 
namical process 1 + X ( A + B ) ; as is easily checked. On the contrary, Trotter formula |TJ] 



exp [x(A + B)] 
based on the identity 

exp [x(A + B)] = lim 



exp( — )exp( 

n n 




(7) 



, xA . . xB ^ 
exp( — ) exp( 



n 



a 



preserves the unitary property in each elementary process exp(^) exp(^-), as is also 
easily checked. The second order formula called Leap flog method has the form 



exp [x(A + B)} 



,xA s 



,xB s 



.xA s 



exp( — J exp( — ) exp{ — 
2n n 2n 



+ 



x 



(9) 



Furthermore, many other higher-order formulas for exponential operators exp [x(A + B)] 
preserving the symmetry corresponding to the unitary propery were recently discovered 
independently both in the development of simulation technology called quantum Monte 
Carlo methods |TU| , |T5| ] to simulate density matrices, or in the development of simula- 
tion technology called symplectic integrators |9|, [11], [TJ, [13], p~|, [23|| to simulate classical 
Hamiltonian dynamical systems. It is an easy task to extend these decomposition for- 
mula of the exponential operators exp [x(A + B)] to more generalized exponential oper- 
ators exp xY}j=iAj of multi noncommutative operators A 1 , A 2 , ■ ■ ■ , Ai. Thus, an ap- 
plication of these successive composition formulas of exponential operators to quantum 
computations of exp x Yjj=i Aj can S^ ve us a unified view of this kind of simulations 
as follows: Let us consider the problem of approximately simulating exp x Y^j=\ Aj for 
t < t' < t + At based on an explicit algorithm on quantum model of computation whose 



3 



each elementary process is successively generated by explicitly time-dependent Hamilto- 
nians Qj(t, At), 1 < j < m. Then, each s-th order approximation formula has a form: 



exp 



i 



[I exp [xQj(t,At)] + 0(z 



■s+l\ 



(10) 



3=1 



The relation 



3=1 3=1 



'ID 



must hold from the lowest order terms in x in Eq. (110). 



3 Theorem 

Let us consider a time-independent Hamiltonian H(q,p) in a certain class of the set 
of Hermitian operators G = {G(q,p)}, where q and p denotes the canonical conjugate 
operators in the standard sense of quantum mechanics. We can define quantum non- 
integrability as follows: 



Definition 1 We call a quantum Hamiltonian system with a time-independent Hamilto- 
nian operator H quantum non-integrable if the following relation holds: 

[$, H] = =>- $ = F(H), (12) 

where $gG and F is a some function of a variable. 

Since if is a time-independent Hamiltonian operator, the expectation value of H must 
be preserved: 

4 < H >= 4 < >= 0, (13) 

at dt 

where < \&| is the state vector. 

Here, we prove the following theorem: 



Theorem 1 If an explicit algorithm preserving unitary property can simulate a quan- 
tum non-integrable system with a time-independent Hamiltonian H approximately, the 
conservation law of the expectation value of the Hamiltonian operator < H > must break 
down. 
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Remark 1: This theorem does not depend on the order and types of approximate algo- 
rithms we choose. 

Remark 2: A class of explicit algorithms preserving unitary property involves universal 
quantum Turing machines in the sense of Deutsch[[|. Thus, as is shown in Ref. |19, ^(J 



this theorem shows that there is no (discrete time) quantum computers to simulate quan- 
tum non-integrable systems without breaking the conservation law of the energy expecta- 
tion. However, the present theorem says not only the limitation of quantum computers on 
this aspect but also a more general statement that there is a universal trade-off between 
the unitary property and the conservation law of energy expectation in time-discretization 
of quantum non-integrable systems. 

(Proof of Theorem 1) 

By using the expression of quantum algorithms in Eq. fllQD , we can consider an s-th order 
algorithm of approximately simulating the quantum dynamics of H for the time duration 
At of a computational step as follows: 

m 

exp [xH] = [] ex P [xQj{t, At )} + 0(x s+1 ), (14) 

3=1 

where x = —iAt/h and 1 < s < oo. Each quantum algorithm Qj(t, At) has a correspond- 
ing time- dependent Hamiltonian Hj(t) satisfying 

rt+At °° i r h r tn-i 

Q j (t,At)=T(exp Hj (s)ds) = l + J2(-*r) n dt x --- dt n Hj(t x ) ■ ■■H j (t n ), 

J t ?2=1 

(15) 

where T denotes the time ordering. The resulting quantum algorithm has also an time- 
dependent Hamiltonian H(t, At) satisfying the relation 

m 

JJ exp [xQj(A, t)] = exp(xH(t, At)). (16) 

3=1 

By successively applying the Baker-Campbell-Hausdorff formula: 

expXexpF = expZ, (17) 

where 

Z = X + Y + \[X, Y] + 1 ([X, [X, Y]] + % [Y, X}}) + [Y, [Y, X}}} + ■■■ (18) 

to the system (|T6|), we can compute the corresponding time-dependent Hamiltonian H in 
a form: 

oo 

H{q,p,t,At) =H + Y l (At) a H a (t) = H + 0(x s ), (19) 
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where H s (t) is a time-dependent correction term of order s. We assume that the energy 
expectation < H > in the quantum simulation is also preserved: 



< H >=< H >= Const, for t < t' < t + At. 



(20) 



Since we can choose At an arbitrary real number, the relation (|20"D means the following 
commutation relations hold: 



H, H(t) =0 and [H, H n (t)] = for n > s. 



[21) 



However, from the assumption of quantum non-integrability of H, it follows that that 
H = F(H). This means that the quantum algorithm H generates the exact quantum 
dynamics of H. This exactness (s — > oo) contradicts the assumption that the underlying 
quantum algorithm gives an approximate tracking of the dynamics of H in the finite 
order s. 

(End of Proof) 



4 Discussions 

The key of the present analysis is in quantum non-integrability. How generic is the notion 
of quantum non-integrability in quantum mechanics? In classical mechanics, it is known 
that most dynamical systems are non-integrable since the famous Poincare theorem in the 
last century. Furthermore, we have exact criteria of classical non-integrability for explicitly 
given Hamiltonian systems based on the singularity analysis |§, [16], [L7|, [0|, |22|, |24|, |25[ . 

On the contrary, in quantum mechanics, we do not have any theorem guaranteeing the 
generic character of quantum non-integrable systems corresponding to the Poincare the- 
orem in classical mechanics nor exact criteria of quantum non-integrability for explicitly 
given Hamiltonian operators. In other words, it is not a trivial thing to connect classical 
non-integrability with quantum non-integrability |], fJT|. Recently, the present author 
found that the quantum Hamiltonian system with a time-independent Hamiltonian oper- 
ator H = \{p1+pl + would be quantum- nonintegrable under the hypothesis of the 
Weyl rule for canonical variables Pi, qi using the Moyal bracket, based on Ziglin's result of 
proof of its classical non-integrability |2(| . For the classical system of this system, it was 



shown in Ref. [0 that we cannot avoid energy fluctuations for some specific initial con- 
ditions like (qi,q2,Pi,P2 = (1000,0.002,0,0)) because the higher-order correction terms 
H s also become bigger as 

\H S \ w AB S , (22) 

where A and B are some positive real constants. It can be easily predicted that this diver- 
gence of the higher-order correction terms H s can also occur in quantum non-integrable 
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systems like a quantum version of the above system. This model can be a vivid ex- 
ample causing rather general phenomena of the breakdown of the conservation law of 
energy expectation for quantum non-integrable systems by using any finite-order time- 
descretization preserving unitary property, which Theorem 1 asserts. This result has an 
interesting implication concerning the usual energy-time uncertain relations. From the 
energy-time uncertain relations At • AE > h, it follows that 

At s \H s \ AE > (23) 

This inequality gives a lower-bound of At which depends only on the order s of time- 
discretization and the correction terms H s . Thus, this analysis shows that our naive view 
that the continuous nature of time in quantum mechanics is naturally obtained in the 
continuum limit At — > is not universal, at least in quantum non-integrable systems. 
It will be an interesting and important open problems to consider time-discretization of 
quantum non-integrable systems in connection with the foundation of quantum mechanics. 
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